Quantised orbital angular momentum transfer and magnetic dichroism in the 
interaction of electron vortices with matter 



Sophia Lloyd, Mohamed Babiker, and Jun Yuan 
Deaprtment of Physics, University of York, Heslington, York, YOlO 5DD, UK 

Following the very recent experimental realisation of electron vortices, we consider their interaction 
with matter, in particular the transfer of orbital angular momentum in the context of electron energy 
loss spectroscopy, and the recently observed dichroism in thin film magnetised iron samples. We 
show here that orbital angular momentum exchange does indeed occur between electron vortices 
and the internal electronic-type motion, as well as center of mass motion of atoms in the electric 
dipole approximation. This contrasts with the case of optical vortices where such transfer only 
occurs in transitions involving multipoles higher than the dipole. The physical basis of the observed 
dichroism is explained. 

PACS numbers: 41.85-p, 42.50Tx 



The suggestion for the existence of electron vortices 
(EVs) was first put forward by Bliokh et al. [T]. In 
analogy with optical vortices (OVs) [2HS], EVs are also 
endowed with the property of quantized orbital angu- 
lar momentum (0AM) of ITi per electron, where I is the 
winding number, but differ from OVs in that EVs have 
electric charge, mass, and spin of i. Bliokh et a/.'s highly 
significant prediction was followed by experimental work 
which succeeded in the generation of EVs, beginning with 
the first experiments by Uchida and Tonomura [B^, who 
used a stepped spiral phase plate, followed by the work of 
Verbeeck et al. [7j, who used a binary holographic grat- 
ing with a Y-like point defect. More recently, using this 
holographic method, McMorran et al. have shown that it 
is possible to generate beams with winding numbers as 
high as Z = 100 [8 . It is now clear that EVs are easily 
generated in conventional electron microscopes with the 
help of holographic reproduction techniques. The poten- 
tial applications are both diverse and wide-ranging and 
it is now known that EVs are set to revolutionize elec- 
tron microscopy and spectroscopy, promising markedly 
improved sub- nanometre resolution in the imaging of ma- 
terials, including biological specimens of materials with 
low absorption contrast [S]. 

As shown theoretically by Berry [TD] and demonstrated 
amply in the context of OVs [2]-[S], it is clearly reason- 
able to state that the orbital angular momentum is a well 
defined quantised property of EVs and hence should be 
exchanged when an electron vortex beam interacts with 
an atom, a molecule, or a solid; for instance, in electron 
energy loss spectroscopy (EELS) using electron vortex 
beams. Furthermore, one expects OAM to be exchanged 
not just in electric dipole interactions, but also in elec- 
tric quadrupole and higher electric multipole interactions 
involving the bound state internal 'electronic-type' mo- 
tion, as well as the gross 'center of mass-type' motion. 
However, in an electric dipole transition the exchange of 
OAM between OVs and matter has been shown |TT] to af- 
fect only the center of mass-type motion, so that to this 
leading multipolar order no OAM transfer occurs from 



OVs to the internal 'electronic' motion. This prediction 
has been confirmed experimentally [12[ I13j . 

This article reports the results of investigations seeking 
to explore the exchange of OAM between EVs and matter 
and to highlight the respects in which EVs differ from 
OVs in the processes of exchange of OAM in this context. 
We also seek to explain the results of recent experimental 
studies by Verbeeck et al., in which the observation of 
dichroism is reported in electron energy loss spectroscopy 
of thin film magnetized iron samples using EVs [7|. 

For simplicity we concentrate on the most basic model 
for matter, displaying both internal 'electronic-type' 
bound states relative to the center of mass, as well as 
gross motion of the center of mass. The model in ques- 
tion is a two-particle hydrogenic atom interacting with 
the electron vortex beam. The EV beam is also taken 
in the simplest form, namely a Bessel beam, which is 
one of the possible solutions of the paraxial Schrodinger 
equation, having well defined orbital angular momentum 
about the beam axis. The interaction between the beam 
electron and the atom is taken as the leading Coulomb 
interaction. For the total Hamiltonian, we write 

H = H^ + H^ + Hi; + Hinu (1) 

where Hq is the Hamiltonian operator for the unper- 
turbed internal electron motion relative to the center of 
mass; its eigenstates are written as ji/'e) = |n;£;m); the 
well known hydrogenic states such that 

H^\n;e;m)=E„\n;e;m), (2) 

where the symbols have their usual meanings. The sec- 
ond term is the Hamiltonian of the unperturbed 
center of mass motion, with eigenstates written \4'n) = 
|K;L) describing both translational and rotational mo- 
tions of the center of mass; 

H^\K;L) = {Ek + Eb)\K;L), (3) 

where L stands for angular momentum of the center of 
mass about the beam axis; Ej^ and -E^ are the energies 
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associated with the translational and rotational motion 
respectively. The eigenfunctions of the center of mass are 
products of plane waves of wavevector K - proportional to 
exp (iK • R) where R is the center of mass position vari- 
able - and rotational states proportional to exp (iL$/j) 
where is the azimuthal angle of the center of mass in 
cylindrical coordinates relative to the center of the elec- 
tron vortex beam. For the unperturbed EV, the Hq has 
eigenfunctions written as {tPb) = \kz'ikp;l) such that 

\k,; kp- 1) ^ {Ek + Ep) %-kp- 1) , (4) 

where the total energy consists of both translational ax- 
ial motion of axial wavevector kz and rotational motion 
about the beam axis. The explicit eigenstates are given 
in cylindrical coordinates in the form 

where Ai is a normalisation factor and kp is the radial 
component of the wavevector, such that total wavevector 
k"^ — kl + hp. Finally Hint is the interaction Hamilto- 
nian which, to leading order, is given by the Coulomb 
interaction energy 

where Cq = e^/(47reo) and qe defines the coordinates of 
the atomic electron relative to the center of the beam. 
In writing Eq.Q, we have made the usual simplification 
that the nucleus is much more massive than the electron 
so that the center of mass is essentially located at the 
nuclear coordinate, R. Note that, in the above, the az- 
imuthal angles denoted by the capital $ are measured 
with respect to the cylindrical coordinate frame centered 
about the electron beam. Azimuthal angles denoted by 
(j) are measured with respect to the spherical coordinate 
frame positioned at the nuclear center of mass. 

The evaluation of effects arising from the interaction 
between the EV state (henceforth to be referred to as the 
beam electron) with the two-particle atom requires con- 
sideration of the transition matrix element Mif where 
|z) = IV'e; V'w; V's) and I/) = IV-e; V'at; V's)- We have 

■Mif = {■llj'^;ip'pf;iljB\Hznt\lpe;ipN;tpB) ■ (7) 

An important technical problem we have to consider here 
is that the atomic electron wavefunctions are normally 
given in spherical polar coordinates relative to the atomic 
center of mass R, whilst in Eq.([6]) both the center of 
mass and the atomic electron are expressed relative to 
the center of the EV beam, using cylindrical polar co- 
ordinates. We are interested in evaluating the transi- 
tion matrix element in such a manner that information 
can be gained as to whether, due to the interaction, a 
transfer of orbital angular momentum has occurred be- 
tween the three subsystems; namely the EV beam, the 



atomic 




FIG. 1. The relevant coordinate frames in the description of 
the interaction between a Bessel electron vortex and a two- 
particle neutral system. Position within the beam, r, is de- 
scribed in cylindrical coordinates about the frame xyz; the 
center of mass of the atom is located at R, and the atomic 
electron is descibed by q in spherical coordinates relative to 
the center of mass. Projection in the xy plane indicates az- 
imuthal angles and of the beam, atomic center of 
mass and atomic electron respectively. 

atomic center of mass and the atomic electron. Since 
the transitions between internal atomic states are classi- 
fied as electric dipole-active, electric quadrupole-active, 
etc., we seek to determine the transfer selection rules in 
the different multipolar approximations. We wish to find 
out whether an experiment can be performed to detect 
0AM exchange between the electron vortex beam and 
the internal dynamics of an atomic or molecular system 
through changes involving electric dipole transitions. It 
turns out, and this constitutes one of the main findings 
of this Letter, that in contrast to the case of optical vor- 
tex interaction with such systems, as discussed in [TT], 
the transfer of 0AM from electron vortices is possible in 
the electric dipole approximation. We shall therefore re- 
strict ourselves to this approximation in the interaction 
Hamiltonian, Eq.([6]). We have 



where q = (qe — R) is the internal electron coordinate 
relative to the center of mass, such that eq is the electric 
dipole moment. 

To proceed we compute explicitly the scalar product in 
Eq.(|8]). With reference to Figjl] and after simplification. 
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the interaction Hamiltonian can be written as 



Hint 



^2 IXr COs((/)e - ^r) - QXR COs(0e - ^ B.) + 

[J--gcos($,-$fl)]5 



(9) 

where Xr, XRj , and Q are functions independent 
of the electron-center of mass coordinate q. In the next 
step towards evaluating the matrix element Eq.([7]), we 
evaluate the integral over azimuthal angles. We write 



X e^^^-^'^-^^iJ^td^.d^H^e, (10) 

where m and m! are magnetic quantum numbers associ- 
ated with the transition between atomic states |ri; £; m > 
and m' >. In evaluating this integral we make the 

substitution j/ = $r ^ ^r, in order to isolate the integral 
over We find that this integral can be expressed in 
terms of generic integrals of the form 



27r ^i(l-l'+a)v 

[J--ecos(y)]3/2 



(11) 



where a is an integer. The final result emerging from the 
angular integral can be written as 

A^azimuth OC Q (5[(;+L) ,(('+L' + l)] <5m,m' -1 

+ '5 ^l(l+L).(l' + L'-l)]^m,m' + l 

+ 1^ 5[(l+L),(l' +L')]5m,m' , (12) 

where Q, S and lA are complex functions of coordinates 
and are all first order in the internal atomic electron co- 
ordinate q, consistent with electric dipole transitions. In 
particular, it turns out that 



Q* =S. 



(13) 



The interpretation of the result in Eq.(12| is as follows 



The first term indicates a dipole active transition in 
which the EV beam or the atomic center of mass loses one 
unit of orbital angular momentum due to a upwards tran- 
sition in which the atom gains one unit of orbital angular 
momentum. The second term indicates that the one unit 
transfer occurs from the atomic electron to either the EV 
beam or the center of mass. The final term shows that 
transitions can occur in which transfer of orbital angular 
momentum occurs only between the vortex and the cen- 
ter of mass rotation, and the internal electronic motion 
is not involved. In situations where the atomic center 
of mass is rigidly fixed, as in a solid, so that L — L' , 
the exchange of orbital angular momentum would occur 
between the EV and the atomic electron. Thus we have 
shown that the electric dipole transition can mediate the 
transfer of a single unit of orbital angular momentum be- 
tween the vortex and the internal motion of the atomic 
electron. Similarly we can show that higher multipole 



transitions lead to exchange of two or more units of or- 
bital angular momentum between the vortex beam and 
the atomic system, but this will not be further discussed 
here. 

The difference in dipole selection rules between OVs 
and EVs - that absorption is forbidden for an OV and 
allowed by an EV - is interesting, as it suggests that an 
0AM electron beam could be more useful in the con- 
text of spectroscopy, as the dipole transition is often the 
dominant process in most physical systems. For exam- 
ple, it means that EV beams could be used to detect 
circular dichroic activity in biological molecules such as 
proteins, in order to gain information about their sec- 
ondary structures [HI [15]. If radiation damage effects 
can be mitigated, the sensitivity and spatial resolution 
will be high. 

The selection rules obtained above form the basis for 
the explanation of the origin of the dichroic signal ob- 
served in the L2 and L3 magnetized iron thin film by 
Verbeeck et al. [7]. As only the m-selection rule is of 
importance in this single particle transition matrix el- 
ement, the same selection rules will emerge if we re- 
place the single particle wavefunction of the hydrogenic 
atom in the above example by the many-particle wave- 
functions of a transition metal atom, i.e. replacing 
\i) = \n-l-m) by |2p63d"; j; m,), and |/) = \n'-l'-m') by 
\2p^2><F'^^ ; j' ]m'j) , where j and rrij are the orbital angu- 
lar momentum and magnetic quantum numbers charac- 
terizing the many-particle configuration consisting of the 
participating 2p core electrons and id valence electrons 
in the transition metal atom [16' . 

Using the matrix elements defined above we first con- 
sider the azimuthal angular dependence of the transition 
rate Ti for transition induced by an EV having / = ±1, 
using Fermi's golden rule. 



F; OC iTWazimuthPp/, 



(14) 



where p/ is the density of final states. Once again only 
the angular factor is needed for our purpose, as we are 
interested in the effect of the beam's orbital angular mo- 
mentum on the transition. The relevant transitions are 
the to' = m ± 1 transitions, and in order for these tran- 
sitions to occur, the one unit of angular momentum of 
I — ±.1 must be gained or lost by the electron vortex. 
We therefore seek to compare the transition rates for I 
and ~L Consider first the dipole transition m' = to + 1, 
I = +1. We find from Eq.(IT2| 



M 

while for to' = m — 1. I 



i=+i 

azimuth 



cx Q, 



T we have from 



"'^azimuth ^ 



-s 



-Q*, 



(15) 



(16) 



where the same proportionality factor applies in the two 
cases. This proportionality factor includes the radial in- 
tegrals, involving Bessel functions arising from the initial 
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and final beam states, and since J-i{x) = — J-|-i(a;) a 
factor of —1 arises in Eq.(16|. It thus follows, since the 
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transition rate depends on the modulus square of the 
matrix element which is the same for / = ±1, and pro- 
vided that the densities of final magnetic states remain 
the same, i.e. Pm+i = Pm~i, we will have 



r, 



(17) 



This means there is no difference in the strength between 
the two transitions involving different helicities I — ±1. 
There is therefore no dichroism for non-magnetic transi- 
tions. 

The clear dichroism observed in the experiment by Ver- 
beeck et al. can thus be explained as arising due to 
the difference in the density of states p/, which will differ 
for each m! . The fact that it is a dipole mediated transi- 
tion means that the dichroic signal can be measured us- 
ing small angle scattering, where such dipole interactions 
dominate. This has clear advantages compared to the 
electron energy-loss magnetic chiral dichroism (EMCD) 
technique devised by Schattschneider et al. [T7| HH] in 
which large scattering angles are required. 

In conclusion, we have shown by direct analysis that it 
is possible to transfer 0AM between an electron vortex 
beam and the internal electron states of an atom in the 
dipole transition and we have checked by direct analy- 
sis that 0AM transfer occurs for quadupole transitions 
and in principle in the case of all higher multipoles. This 
is in direct contrast to the case of optical 0AM trans- 
fer in the interaction with similar systems. It has been 
demonstrated both theoretically and experimentally that 
optical vortices are not specific in their interaction with 
chiral matter. Here we have shown that although orbital 
angular momentum transfer can occur between electron 
vortices and matter in electric dipole transitions for a 
given topological charge |/|, there is no intrinsic differ- 
ence in absorption for the two opposite helicities ±/. We 
have concluded that the 0AM dichroic electron energy- 
loss spectroscopy of the type performed by Verbeeck et 
al. [7] shows a dichroism due to the magnetic nature of 
the material, in which magnetic sublevels would be un- 
equally populated. 
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